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(57) Abstract: The invention relates 
to a method of processing a sound 
field representation which can be used 
to provide a representation which 
has been modified by means of at 
least spatial processing. The inventive 
method comprises the following 
steps: the establishment (8) of an 
initial representation of the sound 
field and the establishment (10) of 
a representation which has been 
modified by filtering combinations 
corresponding to the aforementioned 
processing operation. The invention is 
characterised in that it also comprises 
a step (2) which consists in defining 
the processing operation and which 
comprises: a step (4) involving the 
establishment of a set of directivity 
functions, which is representative of the 
aforementioned processing and which 
takes the form of a set of coefficients 
corresponding to the spherical harmonic decomposition of said functions; and a step (6) in which the above-mentioned filtering 
combinations are determined from the decomposition so that another step (10) involving the application of the thus defined filtering 
combinations can be performed. 

(57) Abrege : Ce procecle" de traitement d'une representation d'un champ acoustique permet de deiivrer une representation modiflee 
par un traitement au moins spatial et comporte :- retablissement (8) d'une representation initiale du champ acoustique; et- 1'eta- 
blissement (10) d'une representation modiflee par des combinaisons de filtrage correspondant audit traitement. H est caracterise* 
en ce qu'il comporte en outre une 6tape (2) de definition dudit traitement, laquelle comprend :- une etape (4) d'etablissement d'un 
ensemble de fonctions de directivite representatif dudit traitement, sous la forme d'un ensemble de coefficients correspondant a la 
decomposition de ces fonctions sur une base d'harmoniques spheriques ; et - une etape (6) de determination desdites combinaisons 
de filtrage a partir de cette decomposition pour la mise en uvre d'une etape (10) d' application des combinaisons de filtrage ainsi 
definies. 
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Procede et systeme de traitement d'une representation 
d 9 un champ acoustique. 

La presente invention concerne un proc6d6 et un systeme de traite- 
ment d'une representation d'un champ acoustique afin de deiivrer une represen- 
tation de ce champ acoustique modifiee par un traitement au moins spatial. 

Actuellement, il existe des proced6s et des systemes qui, & partir 
5 d'une representation d'un champ acoustique sous la forme de coefficients en 
fonction du temps et de I'espace, permettent de deiivrer une representation modi- 
fiee de ces coefficients par des combinaisons lineaires correspondant d un trai- 
tement au moins spatial du champ acoustique. 

Cependant, les traitements possibles sont limites et dans tous les cas 
1 0 predetermines. 

Par exemple, les systemes regroup6s sous le terme « ambisonic » uti- 
lisent une representation d'un environnement sonore sous la forme d'une 
distribution angulaire de sources virtuelles autour d'un point assimiie & la position 
d'ecoute. 

15 Ces systemes ne permettent d'effectuer que certaines modifications 

pr6determin6es telles qu'une distorsion angulaire sur I'azimut, une focalisation 
dans une direction privilegiee, ou une rotation par rapport a un axe donne. 

De plus, ces systemes sont limites au niveau de leur precision. Ainsi, 
dans les systemes existants, des traitements tels que la focalisation dans une 

20 direction quelconque ne peuvent etre appliques qu'£ des representations limitees 
a un niveau de precision correspondant au premier ordre mathematique, c'est a 
dire de faible qualite. 

II apparatt done qu'aucun des systemes de I'art anterieur ne permef 
d'effectuer une modification quelconque d'une representation d'un champ acous- 

25 tique afin d'obtenir une representation modifiee, notamment dans le cas d'une 
representation initiale de qualite, d'ordre 6lev6 par exemple. 

Le but de I'invention est de rem6dier d ce probleme en fournissant un 
procede et un systeme de traitement de la representation d'un champ acoustique 
permettant d'effectuer des traitements quelconques, pr6sentant une tres large 

30 diversite 

L'invention a pour objet un proc6d6 de traitement d'une representation 
d'un champ acoustique, afin de deiivrer une representation modifiee par applica- 
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tion de ce traitement correspondant & un traitement au moins spatial du champ 
acoustique, comportant : 

- retablissement d'une representation initiale du champ acoustique 
sous la forme d'un ensemble de coefficients repr6sentatifs de ce champ dans le 

5 temps et dans les trois dimensions de I'espace ; et 

- retablissement d'une representation modifiee par des combinaisons 
de filtrage correspondant audit traitement, appliqu6es auxdits coefficients de la- 
dite representation initiale, 

caracterise en ce qu'il comporte en outre une etape de definition dudit 
10 traitement, laquelle comprend : 

- une etape d'etablissement d'un ensemble de fonctions de directivity 
repr6sentatif dudit traitement, sous la forme d'un ensemble de coefficients 
correspondant a la decomposition de ces fonctions sur une base d'harmoniques 
spheriques ; et 

15 - une etape de determination desdites combinaisons de filtrage a partir 

de cette decomposition pour la mise en oeuvre d'une etape duplications des 
combinaisons de filtrage ainsi definies auxdits coefficients de la representation 
initiale, afin d'obtenir ladite representation modifiee. 

Suivant d'autres caracteristiques de I'invention : 

20 - ladite etape de determination d'un ensemble de fonctions de directi- 

vite comporte une sous-etape de determination directe de tout ou partie dudit 
ensemble de fonctions de directivite dans Tun ou I'autre d'une pluralite de formats 
predetermines ; 

- ladite etape de determination d'un ensemble de fonctions de directivi- 
25 te, comporte une sous-6tape d'interpretation automatique d'une consigne de trai- 
tement permettant de delivrer automatiquement tout ou partie dudit ensemble de 
fonctions de directivite dans I'un ou I'autre d'une pluralite de formats predetermi- 
nes, a partir de ladite consigne de traitement ; 

- ledit ensemble de fonctions de directivite representatif du traitement 
30 est exprime selon Tun ou I'autre d'une pluralite de formats predetermines et ladite 

etape d'etablissement d'un ensemble de fonctions de directivite comporte une 
sous-etape de conversion, a partir dudit format predetermine, en un ensemble de 
coefficients correspondant & la decomposition desdites fonctions de directivite sur 
une base d'harmoniques spheriques ; 
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- ladite etape de definition du traitement comprend en outre, une etape 
d'etablissement d'un ensemble de parametres dudit traitement ; 

- ledit ensemble de parametres comprend tout ou partie des parame- 
tres suivants : 

5 - un ensemble de constantes formant des coefficients de pond6ra- 

tion et defini en fonction de la nature dudit traitement ; et 

- un parametre specifiant I'ordre limite de fonctionnement du prece- 
de et correspondant a la limite mathematique de precision souhaitee ; 

- ladite etape d'etablissement d'un ensemble de parametres comporte 
10 une sous-6tape de determination directe de tout ou partie dudit ensemble de pa- 
rametres ; 

- ladite etape d'etablissement d'un ensemble de parametres comporte 
une sous-6tape d'interpretation automatique d'une consigne de traitement per- 
mettant de deiivrer automatiquement tout ou partie dudit ensemble de parame- 

15 tres £ partir de ladite consigne de traitement ; 

- ladite etape d'etablissement d'un ensemble de parametres comprend 
une sous-6tape de calcul de tout ou partie dudit ensemble de constantes en fonc- 
tion de la nature dudit traitement ; 

- ladite etape de determination desdites combinaisons de filtrage com- 
20 prend au moins une sous-6tape de calcul de combinaisons de filtrage par des 

combinaisons Iin6aires desdits coefficients correspondant a la decomposition sur 
une base d'harmoniques spheriques dudit ensemble de fonctions de directivity 
repr6sentatif dudit traitement ; 

- ladite sous-6tape de calcul par des combinaisons Iin6aires fait inter- 
25 venir ledit ensemble de constantes en tant que coefficients de pond6ration ; 

- il est r6pet6 une pluralite de fois et les combinaisons de filtrage resul- 
tant des differentes sous-etapes de calcul sont combin6es entre elles pour former 
les combinaisons de filtrage destinies a etre appliqu6es auxdits coefficients de la 
representation initiale ; 

30 - le resultat de ladite sous-£tape de calcul forme directement lesdites 

combinaisons de filtrage destin6es d §tre appliqu6es auxdits coefficients de la 
representation initiale ; 

- ledit traitement est un traitement de convolution ; 

- ledit traitement est un traitement de fenetrage ; 
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- ledit traitement comprend au moins une distorsion £ appliquer a la 
representation initiale du champ acoustique et tout ou partie des fonctions de di- 
rectivity determin6es lors de retape d'etablissement d'un ensemble de fonctions 
de directivity torment un ensemble de couples de distorsion representatif de la- 

5 dite distorsion sous la forme d'un ensemble de coefficients correspondant d la 
decomposition desdits couples de distorsion sur une base d'harmoniques spheri- 
ques ; 

- ledit traitement comprend un premier traitement dit nominal et la ou 
lesdites distortions, ladite sous-6tape de calcul deiivrant des combinaisons de 

10 filtrage correspondant audit traitement nominal et ladite etape de determination 
de combinaisons de filtrage comprenant en outre au moins une sous-etape de 
distorsion par modification des combinaisons de filtrage deiivrees par ladite sous- 
etape de calcul, a partir desdits couples de distorsion repr&sentatifs de la ou des- 
dites distorsions, afin de delivrer lesdites combinaisons de filtrage destinies a 

15 etre appliqu6es auxdits coefficients de la representation initiale ; 

- ledit ensemble de parametres comporte un parametre representatif de I'effort de 
distorsion, soit encore sp6cifiant la ressemblance entre ledit traitement d appli- 
quer et ledit traitement nominal ; 

- ledit traitement comprend uniquement la ou lesdites distorsions, 

20 ladite etape de determination des combinaisons de filtrage comprenant au moins 
une sous-etape de distorsion deiivrant des combinaisons de filtrage uniquement 
a partir desdits couples de distorsion repr6sentatifs de la ou desdites distorsions, 
afin de delivrer lesdites combinaisons de filtrage destinees a etre appliquees 
auxdits coefficients de la representation initiale ; 

25 - il comporte une plurality de sous-etapes de distorsion, r6cursives les 

unes par rapport aux autres, de sorte que chaque sous-etape de distorsion regoit 
en entree les combinaisons de filtrage delivrees par la sous-etape precedente ; 

- ladite au moins une sous-6tape de distorsion comprend une sous- 
etape d'optimisation permettant de determiner les combinaisons de filtrage 

30 correspondant d des modifications a effectuer sur des parties de la representa- 
tion initiale qui ne sont pas directement visees par lesdits couples de distorsion ; 

- ledit traitement se decompose en une pluralite de traitements 6I6- 
mentaires et au moins ladite etape de determination de combinaisons de filtrage 
est repetee une pluralite de fois, les combinaisons de filtrage resultant des diffe- 
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rentes sous-etapes de calcul et/ou de distorsion et correspondent auxdits traite- 
ments elementaires etant arrang6es entre elles pour former les combinaisons de 
filtrage destinees a etre appliqu6es auxdits coefficients de la representation ini- 
tiate ; 

5 - ladite etape de determination des combinaisons de filtrage comprend 

une sous-etape de conversion desdites combinaisons de filtrage representees 
dans un format predetermine pour deiivrer des parametres de filtres adaptes pour 
le parametrage de filtres donn6s en vue de la realisation de tadite etape 
d'application des combinaisons de filtrage ; 
10 - (edit traitement est un traitement spatial et ladite etape duplications 

des combinaisons de filtrage comprend des modifications de gain appliqu6es 
auxdits coefficients de la representation initiate ; 

- ledit traitement est un traitement spatial et fr6quentiel et ladite etape 
d'application des combinaisons de filtrage comprend Implication auxdits coeffi- 

15 cients de la representation initiate, de filtres variant en fonction de la frequence ; 

- lesdits ensembles de coefficients repr6sentatifs de la representation 
initiale du champ acoustique et de la representation modifiee sont des ensembles 
de coefficients dits de Fourier-Bessel. 

[.'invention a egalement pour objet un dispositif de traitement d'une re- 
20 presentation d'un champ acoustique afin de deiivrer une representation de ce 
champ acoustique par application de ce traitement correspondant dun traitement 
au moins spatial du champ acoustique, ledit dispositif recevant en entree une 
representation initiale du champ acoustique sous la forme d'un ensemble de 
coefficients representatif de ce champ dans le temps et dans les trois dimensions 
25 de I'espace et deiivrant en sortie une representation modifiee par les combinai- 
sons de filtrage correspondant audit traitement applique auxdits coefficients de la 
representation initiale, caracterise en ce qu'il comporte des moyens de definition 
dudit traitement, comprenant : 

- des moyens d'etablissement d'un ensemble de fonctions de direc- 
30 tivite representatif dudit traitement, sous la forme d'un ensemble de coefficients 

correspondant & la decomposition de ces fonctions de directivite sur une base 
d'harmoniques sph6riques ; et 

- des moyens de determination desdites combinaisons de filtrage a 
partir de cette decomposition pour leur mise en oeuvre par des moyens (80) de 
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filtrage recevant en entree lesdits coefficients de la representation initiate et deli- 
vrant en sortie ladite representation modifiee. 

Selon d'autres caracteristique de ce dispositif : 

- lesdits moyens de definition dudit traitement comportent un module 
5 de saisie directe d'informations de fonctionnement du dispositif dans Tun ou 

I'autre d'une pluralite de formats ; 

- lesdits moyens de definition dudit traitement comportent un module 
d' interpretation automatique d'une consigne de traitement permettant de delivrer 
automatiquement des informations de fonctionnement dans Tun ou I'autre d'une 

10 pluralite de formats predetermines ; 

- lesdites informations de fonctionnement comprennent ledit ensemble 
de fonctions de directivite repr£sentatif dudit traitement : 

- lesdits moyens de definition dudit traitement comportent un module 
de conversion adapte pour convertir ledit ensemble de fonctions de directivite de 

15 Tun ou I'autre d'une pluralite de formats predetermines en un ensemble de coeffi- 
cients correspondant & la decomposition desdites fonctions de directivite sur une 
base d'harmoniques spheriques ; 

- lesdites informations de fonctionnement comprennent un ensemble 
de parametres dudit traitement ; 

20 - lesdits moyens de definition dudit traitement comportent un module 

de calcul de tout ou partie dudit ensemble de parametres en fonction de la nature 
dudit traitement : 

- lesdits moyens de determination desdites combinaisons de filtrage 
comprennent un module de calcul de combinaisons de filtrage par combinaisons 

25 Iin6aires desdits coefficients correspondant & la decomposition sur une base 
d'harmoniques spheriques dudit ensemble de directivite representatif dudit trai- 
tement ; 

- ledit ensemble des parametres dudit traitement comprend un ensem- 
ble de constantes et ledit module de calcul par combinaisons Iin6aires fait inter- 

30 venir ledit ensemble de constantes en tant que coefficients de pond6ration ; 

- ledit traitement comprend au moins une distorsion, tout ou partie du- 
dit ensemble de fonctions de directivite formant un ensemble de couples de dis- 
torsion correspondant a la ou aux distorsions d appliquer & la representation ini- 
tiate et lesdits moyens de determination desdites combinaisons de filtrage com- 
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prennent en outre un module de distorsion par modification des combinaisons de 
filtrage deiivrees par ledit module de calcul, a partir desdits couples de distorsion 
afin de deiivrer lesdites combinaisons de filtrage destinies a etre appliqu6es 
auxdits coefficients de la representation initiate ; 
5 - ledit traitement comprend uniquement une ou plusieurs distorsions, 

tout ou partie dudit ensemble de fonctions de directivite formant un ensemble de 
couples de distorsion correspondent a la ou aux distorsions & appliquer d la re- 
presentation initiate et lesdits moyens de determination desdites combinaisons de 
filtrage comprennent uniquement un module de distorsion pour deiivrer lesdites 
10 combinaisons de filtrage uniquement & partir desdits couples de distorsion. 

^invention sera mieux comprise a la lecture de la description qui va 
suivre donn6e uniquement a titre d'exemple, et faite en se referent aux dessins 
annexes sur lesquels : 

- la Fig.1 est une representation d'un repere sph6rique ; 

15 - la Fig.2 est un organigramme general du precede de invention ; 

- la Fig.3 est un organigramme detailie d'un premier mode de realisa- 
tion de Tinvention ; 

- la Fig.4 est un schema d'un mode de realisation de retape 
d'application du traitement ; 

20 - la Fig.5 est un organigramme detailie d'un second mode de realisa- 

tion de Tinvention ; et 

- la Fig.6 est un schema synoptique d'un dispositif adapte pour la mise 
en oeuvre du precede de Tinvention. 

Sur la figure 1, on a repr6sent6 un repere sph6rique classique, de 
25 maniere d pr6ciser le systeme de coordonn6es auquel it est fait reference 
dans le texte. 

Ce repere est un repere orthonormal, d'origine 0, et comporte trois 
axes (OX), (OY) et (OZ). 

Dans ce repere, une position notee x est decrite au moyen de ses 
30 coordonn6es spheriques (r 9 0,0), ou r d6signe la distance par rapport a Torigine O t 
^Torientation dans le plan vertical et ^'orientation dans le plan horizontal. 

Dans un tel repere, un champ acoustique est connu si Ton definit en 
tout point d chaque instant t la pression acoustique not6e p{r,0,<f>,t\ dont la trans- 
ferase de Fourier est not6e P(r>0 t <f>j) ou / designe la frequence. 
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Le proc6d6 de invention se base sur Tutilisation de fonctions spatio- 
temporelles permettant de d6crire un champ acoustique quelconque dans le 
temps et dans les trois dimensions de I'espace. 

Dans les modes de realisation d6crits, ces fonctions sont des fonctions 
5 dites de Fourier-Bessel sph6riques de premiere esp6ce, appel6es par la suite 
fonctions de Fourier-Bessel. 

Dans une zone vide de sources et vide d'obstacles, les fonctions de 
Fourier-Bessel correspondent aux solutions de liquation des ondes et consti- 
tuent une base qui engendre tous les champs acoustiques produits par des sour- 
1 0 ces situ6es a I'exterieur de cette zone. 

Tout champ acoustique tridimensionnel peut done s'exprimer par une 
combinaison lineaire des fonctions de Fourier-Bessel, selon I'expression sui- 
vante : 

15 Dans cette equation, les termes /%,(/) sont d6finis comme les coeffi- 

cients de Fourier-Bessel du champ p(r,0,fci), k = ^L, c est la c6l6rit6 du son 

c 

dans Tair (340 ms" 1 ), jikr) est la fonction de Bessel spherique de premiere es- 
p§ce d'ordre / definie par ji(x) = J UV2 (x) ou Jy(x) est la fonction de Bessel de 

premiere esp^ce d'ordre v % et yr{0$ est I'harmonique sph6rique r6elle d'ordre / 
20 et de terme m, avec m allant de -/ a /, d6finie par : 

yr(W= P} m] (cos0)trg m (0) 

avec : 



— t-cos(m^) pour m>0 
trgm(0)=\ } pour m=0 

-4= sin( m <f>) pour m<0 

Dans cette 6quation, les Pr(x) sont les fonctions de Legendre asso- 
25 ctees d6finies par : 

avec Pfe) d6signant les polyndmes de Legendre, d6finis par : 
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Les coefficients de Fourier-Bessel s'expriment aussi dans le domaine 
tempore! par les coefficients pi^i) correspondant a la transform6e de Fourier 
temporelle inverse des coefficients de Fourier-Bessel Pi^tf). 
5 Dans d'autres modes de realisation, le champ acoustique est decom- 

pose sur une base de fonctions, ou chacune des fonctions s'exprime par une 
combinaison Iin6aire, 6ventuellement infinies de fonctions de Fourier-Bessel. 

Sur la figure 2, on a repr&sente un organigramme g6n6ral du proc6de 
de Tinvention. 

10 De maniere gen6rale, le proc6d6 debute par une etape 2 de definition 

d'un traitement & appliquer d une representation d'un champ acoustique et 
correspondant a un traitement au moins spatial de ce champ acoustique. 

Cette etape 2 de definition du traitement comporte une etape 4 
d'etablissement d'un ensemble de fonctions de directivite representatif du traite- 

15 ment. 

On entend par fonction de directivite une fonction, notee Giff^J), qui 
associe une valeur reelle ou complexe, 6ventuellement fonction de la frequence, 
ou d'une plage de frequences, d chaque direction de Tespace. 

Toute fonction de directivite G(0,0J) se decompose en coefficients 
20 d'harmoniques spheriques, notes par transform6e de Fourier spherique. 

Ainsi, la fonction de directivite G{0^J) s'obtient a partir des coefficients de Fourier 
spherique Gijj) par transform6e de Fourier spherique inverse, selon 
I'expression : 

/=0m=-/ 

25 Dans le reste du document, les fonctions de directivite peuvent etre 

designees par leur expression G{0 y 0j) dans I'espace direct ou sous forme de 
coefficients de Fourier spherique G/^,(/). 

Cet ensemble de fonctions de directivite est deiivre a Tissue de retape 
4 sous la forme d'un ensemble de coefficients correspondant a la decomposition 
30 des fonctions de directivite sur une base d'harmoniques spheriques. 

Par la suite, une etape 6 de determination de combinaisons de filtrage 
est realisee a partir de la decomposition des fonctions de directivite realisee & 
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l'6tape 4. La notion de filtrage recouvre le filtrage frequentiel, d'amplitude, de 
phase ou leur combinaison. 

Ainsi, I'etape 6 de determination de combinaison de filtrage deiivre des 
combinaisons de filtrage correspondant au traitement a appliquer sous la forme 
5 de parametres de filtres permettant ulterieurement de mettre en oeuvre ces com- 
binaisons de filtrage a I'aide par exemple d'un ensemble de filtres num6riques. 

En vue de I'application du traitement, le procede comporte une etape 8 
d'etablissement d'une representation initiate du champ acoustique sous la forme 
d'un ensemble de coefficients ou signaux representatif de ce champ dans le 
1 0 temps et dans les trois dimensions de I'espace. 

Notamment, cet ensemble de coefficients est un ensemble de coeffi- 
cients dit de « Fourier-Bessel » tels que d6finis en reference a la figure 1. 

Le proc6d6 comporte ensuite une etape 10 duplication des combi- 
naisons de filtrage, definies lors de retape 2 de definition, aux coefficients de la 
15 representation initiale tels que definis lors de I'etape 8. 

Cette etape 10 permet d'obtenir une representation modifiee corres- 
pondant a une representation dans le temps et dans les trois dimensions de 
I'espace du champ acoustique initial modifi6 par le traitement defini, sous la 
forme d'un ensemble de coefficients dits de Fourier-Bessel. 
20 Sur la figure 3, on a repr6sente en detail Torganigramme d'un premier 

mode de realisation de I'invention. 

Dans ce mode de realisation, les traitements appliques sont des trai- 
tements tels que le fenetrage spatial conslstant d modifier le gain de tout ou par- 
tie des sources virtuelles dont on considere qu'elles generent le champ acousti- 
25 que initial, et la convolution spatiale ou filtrage spatial, consistant & modifier 
retalement angulaire de tout ou partie de ces sources virtuelles, ou tout autre 
traitement analogue. 

Ces traitements de fenetrage et convolution peuvent varier avec la fre- 
quence, ce qui permet d'effectuer des modifications d'6loignement et de spectre 
30 de tout ou partie de ces sources virtuelles. 

Ainsi que cela a 6t6 decrit d la figure 2, le procede debute par une 
6tape 4 de determination d'un ensemble de fonctions de directivite representatif 
du traitement £ appliquer. 
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Dans le cas 0C1 le traitement ne comporte qu'un seul traitement tel 
qu'un traitement de fenetrage ou de convolution, une seule fonction de directivity, 
6ventuellement variant avec la frequence, suffit £ le representee 

Avantageusement, I'etape 4 comporte une sous-etape 12 de 
5 determination directe de tout ou partie de cet ensemble de fonctions de directivite 
dans Tun ou I'autre d'une plurality de formats predetermines. 

En outre ou altemativement, I'etape 4 comporte une sous-etape 14 
d'entr6e d'une consigne de traitement, laquelle est interpretee automatiquement 
lors d'une sous-etape 16 d'interpretation permettant de delivrer automatiquement 
10 tout ou partie dudit ensemble de fonctions de directivite dans Tun ou I'autre d'une 
plurality de formats predetermines. Par exemple, le traitement est mis en ceuvre 
par des moyens logiciels sur un ordinateur et cette consigne de traitement a la 
forme d'une ic6ne dont la validation entraTne I'execution par les moyens logiciels 
d'un traitement donne. 
15 Les sous-6tapes 12 et 14 correspondent par exemple a une saisie 

manuelle ou encore d la reception de donn6es de directivite provenant d'une 
memoire interne ou d'un dispositif exteme. 

En fonction du format dans lequel est determine ledit ensemble de 
fonctions de directivite, retape 4 peut requ6rir une sous-etape 18 de conversion 
20 de I'ensemble de fonctions de directivite £ partir d'un format predetermine en un 
ensemble de coefficients correspondant a la decomposition des fonctions de di- 
rectivite sur une base d'harmoniques sph6riques. 

Pour chaque fonction de directivite Gt^Jf) de I'ensemble de fonctions 
de directivite, cette sous-6tape 18 de conversion peut etre realis6e de plusieurs 
25 manieres. 

Par exemple, la fonction de directivite GiJJ) est definie d partir de pa- 
rametres (r g , 0 q> 0 q ) sp6ciflant la position de Q sources virtuelles et de parametres 
H q {f) sp6cifiant la r6ponse en frequence de chacune de ces sources. A chaque 
source virtuelle est associe un modele de rayonnement de sorte que I'ensemble 
30 de ces sources est suppose produire un champ acoustique dont la representation 
en coefficients de Fourier-Bessel permet d'associer £ ce champ, une fonction de 
directivite. 
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Par exemple, le module de rayonnement associe & chaque source est 
un modele de rayonnement en ondes sphgriques et la fonction de directivity 
Gi^ij) est obtenue lors de la sous-6tape 18 gr§ce aux 6quations suivantes : 

5 avec 



En variante, ce calcul est realise a partir de paramytres (0 g9 fa) sp6ci- 
fiant les directions des Q sources. Dans ce cas, un module de rayonnement en 
ondes planes est associ& & chaque source et la fonction de directivity Gi^JJ) est 
10 obtenue gr§ce £ liquation suivante : 

Dans d'autres modes de realisation, G itm (f) est determine £ partir d'une 
combinaison de modules de rayonnement en ondes planes et en ondes spheri- 
ques 

15 Dans un autre cas, la fonction de directivity est d£finie & partir de pa- 

ramytres repr6sentatifs d'une description paramytrique d'une fonction de directi- 
vity. Cette description s'effectue dans I'espace direct et associe & chaque direc- 
tion une valeur pouvant etre complexe. Avantageusement, les paramytres peu- 
vent varier avec la frequence. 

20 Dans ce cas, les coefficients de la fonction de directivity G lim (j) sont dy- 

terminys par la mise en oeuvre d'un modyie paramytrique. Les paramytres du 
modyie sont placys dans un vecteur j{f) et le modyie paramytrique de la fonction 

de directivity s'ycrit G(0,0 9 p(J)). Les coefficients de la fonction de directivity Gijj) 

s'obtiennent par transformye de Fourier sphyrique de G{0^J{f)) selon 

25 ryquation : 

GW/)= P G(^,7(/)) yj»(*,*) sin0 d6d<f> 

Par exemple, le modyie paramytrique peut correspondre a une fonc- 
tion de directivity prysentant un unique lobe. Ce modele est controle par les pa- 
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rametres (#(/),$(/)) sp6cifiant la direction du lobe et le param6tre 6^{f) specifiant 
Tangle d'ouverture du lobe. 

Ainsi, le modele parametrique peut etre une fonction de directivity va- 
lant 0 pour les directions quelconques et valant 1 pour chaque direction v6ri- 
5 fiant la relation suivante : 

(0iV)>tiV)) .(^^>cos(^(/)) 
ou (0[V),0\(J)) . (0,0) d6signe le produit scalaire entre les directions et 



20 



Dans ce cas les coefficients G/^(/) de la fonction de directivity sont 
10 donn6spar: 

Dans cette expression P. x (x) = 1. 

En variante, les parametres de la description param6trique de la fonc- 
tion de directivite peuvent varier avec le temps. Dans ce cas, les parametres du 
15 modele sont places dans un vecteur j{t) et le modele parametrique de la fonction 

de directivity s'ecrit g(0 9 0 9 j{t)). Les coefficients de la fonction de directivity gi^it) 
exprimes dans le domaine temporel s'obtiennent alors par transformee de Fourier 
sphyrique de g(0 9 0 9 j{t) ) selon I'expression : 



9i,m(t) = / / g(9,<f>tf(t)) y?(e,<t>) sm6d0d<f> 

J 9=0 J 6=0 



Les coefficients G /#m (/) de la fonction de directivity exprimys dans le 
domaine fryquentiel s'obtiennent alors par transformee de Fourier temporelle des 
coefficients gisnit). 

Dans encore un autre cas, la fonction de directivity est dyfinie sous la 
25 forme d'une liste de Q echantillons. La liste est formee de Q couples {((0 q9 # q ), v q )} 
ou (0 q ,0q) est la direction de rychantillon q et od v q est la valeur de la fonction de 
directivity dans la direction (0 q ,0 q )- Avantageusement, la liste d'ychantillons peut 
etre saisie pour une plurality de frequences / et est alors notye {((0 q ,0 q ), v q )}(f). 

Pour chaque fryquence / de fonctionnement, les coefficients Gi^if) de 
30 la fonction de directivity sont dyduits de la liste d'ychantillons {((4>A)> v q )}(f). Ces 
coefficients sont obtenus par inversion du processus d'ychantillonnage angulaire 
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qui permet de d6duire les echantillons de la liste {((^,^), v q )}(f) 3 partir d'une 
fonction de directivity fournie sous forme de coefficients d'harmoniques sph6ri- 
ques. Cette inversion peut prendre diff6rentes formes afin de controler 
('interpolation entre les echantillons. 

Dans le mode de realisation decrit, les echantillons sont places dans 
un vecteur Vde la fagon suivante. 

[v, v 2 — v*]' 

Une matrice Y est constitu6e d partir des directions d'echantillonnage 
de la fagon suivante : 





Y l,0,0 ^5,1 ,- 


1 ^1,1,0 ^1,1.1" 




-L" 


' Y u.,o" 


' Y l,L,L 


10 


^2,0,0^2,1,- 


■I ^2,1,0^2,1,1" " 

• « 


' Y 2,L,- 


-V 


' Y 2 K L,0" 


• Y 2 K L,L 




YqmbXqa.- 




' Y Q,L, 


-L" 


' Y Q,L.0~ 


' Y Q,L,L. 



ou chaque 6l6ment Y qJim est 6gal a : 

On rappelle que yF(0 q ,fc) d6signe I'harmonique spherique d'ordre / et 
15 determem pour la direction d'echantillonnage 0 q ,<fa consideree. 

On determine un vecteur G contenant les coefficients pour la 

frequence / Le vecteur G peut etre determine au moyen des relations suivantes : 

- si le nombre Q d'echantillons a la frequence / est superieur au 
nombre de coefficients G iim {f) : 

20 G = (Y T Y)' 1 Y T V\et 

- si le nombre Q d'echantillons k la frequence / est inferieur au 
nombre de coefficients Gi^if) : 

G=Y T (YY T y l V. 

Dans le vecteur G % les coefficients G tim {f) sont organises de la maniere 

25 suivante : 

[GoM d.-iC/) GM GM - G l -lV) - GM - G LtL {f)\ 
En variante, la liste d'echantillons peut etre saisie pour une plurality 
d'instants / et est alors notee {((0 q ,0 q ) 9 v q )}(t). 

Pour chaque instant t saisi, les coefficients gi^if) de la fonction de di- 
30 rectivite exprimes dans le domaine temporel sont d6duits de la liste d'echantillons 
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V *)H0- Ces coefficients sont obtenus par inversion du processus 
d'6chantillonnage angulaire qui permet de deduire les echantillons de la liste 
{((4»6)f V *)H0 * P artir d'une fonction de directivite fournie sous forme de coeffi- 
cients d'harmoniques spheriques. 
5 Les coefficients Gi^if) de la fonction de directivite exprimes dans le 

domaine frequentiel s'obtiennent alors par transformee de Fourier temporelle des 
coefficients gi^t). 

Bien entendu, dans le cas ou les etapes 12 ou 16 deiivrent 
directement la fonction de directivite sous la forme de coefficients correspondent 
10 a la decomposition de cette fonction sur une base d'harmoniques spheriques, la 
sous-etape 18 de conversion ne la modifie pas. 

Ainsi, d Tissue de I'etape 4, un ensemble de fonctions de directivite re- 
presentatif du traitement a appliquer est deiivre sous la forme d'un ensemble de 
coefficients correspondant a la decomposition de ces fonctions sur une base 
1 5 d'harmoniques spheriques. 

Par ailleurs, I'etape 2 de definition du traitement comporte egalement 
une etape 20 d'etablissement d'un ensemble de parametres du traitement. 

De meme que dans le cas de I'etape 4, cet ensemble de parametres 
peut etre determine directement lors d'une sous-etape 22 de determination di- 
20 recte et/ou a la suite d'une sous-etape 24 d'entr6e d'une consigne de traitement 
et d'une sous-etape 26 d'interpretation de cette consigne de traitement permet- 
tant de delivrer automatiquement tout ou partie de I'ensemble de parametres. 

Les sous-etapes 22 et 24 correspondent par exemple £ une saisie 
manuelle de certains parametres ou encore £ la reception de donn6es de para- 
25 metres a partir d'une memoire interne ou d'un dispositif externe. 

Notamment, lors de I'etape 20 d'etablissement d'un ensemble de pa- 
rametres, tout ou partie des parametres suivants sont definis : 

- un parametre L specifiant I'ordre limite de fonctionnement du precede 
et correspondant £ la limite mathematique de precision souhaitee ; et 
30 - un ensemble de constantes not6es cf t ^ r ^ formant des coefficients de 

ponderation definis en fonction de la nature du traitement. En particulier, ces 
coefficients sont avantageusement appliques aux fonctions de directivite pour 
obtenir les combinaisons de filtrage pr6cit6es. 
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Certains parametres peuvent etre invariant et etre directement integres 
dans le proc6d6. 

Le parametre L permet au precede de s'adapter au niveau de preci- 
sion de la representation initiale et permet notamment au precede, objet de 
rinvention, d'appliquer un traitement a une representation initiale pr6sentant un 
haut niveau de precision. 

Ainsi que cela a 6te decrit precedemment, les constantes cfcjj).^. peu- 
vent etre determines directement lors de retape 20 par un operateur, £ partir 
d'une memoire ou d'un dispositif interne. Par exemple, elles sont toutes calcuiees 
prealablement, stockees dans une memoire et s6lectionnees en fonction de la 
nature du traitement. 

Altemativement, une consigne d6finissant la nature du traitement est 
saisie lors de la sous-etape 24 et la sous-etape 26 correspond au calcul de tout 
ou partie de cet ensemble de constantes en fonction de la nature du traitement et 
a partir de relations decrites ulterieurement dans le texte. 

Bien entendu, les constantes cf;^ r ^ peuvent §tre d6termin6es & I'aide 

d'une combinaison des solutions decrites. 

Les parametres etablis lors de retape 20, ainsi que les coefficients 
correspondants aux fonctions de directivite etablis lors de retape 4, sont utilises 
lors de retape 6 de determination des combinaisons de filtrage. 

Cette etape 6 de determination comporte une sous-etape 30 de calcul 
de combinaisons de filtrage par des combinaisons Iin6aires des coefficients 
correspondants a la decomposition sur une base d'harmoniques sph6riques de 
Tensemble de fonctions de directivite repr6sentatif du traitement. 

Lors de cette premiere sous-etape 30 de calcul, des combinaisons de 
filtrage, notees C\£{f) , sont d6termin6es a partir des fonctions de directivite G^if) 

et des constantes c££- m . determines en fonction de la nature du traitement et 

formant des coefficients de pond6ration des combinaisons Iin6aires. Les combi- 
naisons de filtrage sont alors determinees de la fagon suivante : 



Dans la relation pr6c6dente, / et m d6signent respectivement Tordre et 
le terme des coefficients de la representation initiale, /* m' designent Tordre et le 
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terme des coefficients de la repr6sentation modifi6e et /", m" designent I'ordre et 
le terme des coefficients de la fonction de directivity. 

Dans le cas ou le traitement est un traitement de type fen&rage, les 
constantes c££. tW . sont d6termin6s par les relations : 



c l '> m ' - I 

C f,m,i",m" ~ * 



1 / J'.m' (1) r,m' (1) x 0 



J',m' (1) 



l',m' (1) 
C lflJL» t m" 



pour m = 0 



i / i',m' (1) _ I',m' (1) x 0 



(1) 



avec 



JL. (J** (2) , J'.m' (2) \ // > q 



/,m' (2) 



pour m" — 0 



,i',m' (2) 



avec 



1 (J',m' (3) r,-m- (3A m' > 0 

71 l c J,m,J",m" + c J,m,t",m'V P OUr m ^ U 



(',-m' (3) 



J'.m' (2) _ 
c l,m,l",m" ~ 



i',0 (3) 
c i,m,J",m" 



pour m' = 0 



1 / i',-m' (3) _ t>,m> (3) \ , Q 

J7f l c J,m,t",m" c (,m,/'',m"; P our m ^ u 



10 avec 



l',m'(3) _rm'-m 1 V(» + 1)(2T + l)(2i*+l) 
Cl,m,l«,m« - V 2 f + V + V , + ! 



pf pi' p<" 
°H-|m| ^fHjmM ^<"+|m"| 

r M pK| p Ki 



2v^F 

I i' I" rfc p*-|m| rk > r k'-\m'\ pjfc" pfc"-|m"| 



et 



E E E (-D 

fc=|m| Jfc'=|m'| k"=\m"\ 



1 pour x = y 
0 pour x ^ y 



et 



et 



n = 



|m| + \m'\ + \m"\ 



W p!(n-p)! 



Ce traitement correspond £ I'application d'un fenetrage spatial avec Til— 
trage temporel. En effet, les coefficients d'harmoniques sph6riques Gijj) corres- 
pondent d la fonction de directivite G(0 9 #f). De meme, les coefficients de Fourier- 
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Bessel P^iJJ) du champ acoustique initial correspondent a des coefficients 
d'harmoniques sph6riques d'une fonction de directivite I*\0M- Ainsi, le fene- 
trage spatial avec filtrage temporel se traduit, pour chaque fr6quence f, par un 
fenetrage de la fonction de directivite F^iQ&f) par la fonction de directivite 
G(0 y 0j), ce qui correspond a une multiplication pour chaque direction de la 
valeur de ces deux fonctions de directivite, les valeurs de ces deux fonctions de 
directivite pouvant etre complexes. 

Ainsi que cela a 6t6 d6crit pr6c6demment, ce traitement applique au 
champ acoustique a pour effet de modifier ramplitude, I'eloignement et le spectre 
de tout ou partie de I'ensemble des sources virtuelles dont on considere qu'elles 
genferent le champ acoustique. 

La fonction de directivite G{0,0j) peut etre ind6pendante de la fre- 
quence et est alors note G(0 9 0). Le traitement est alors uniquement spatial et 
correspond, pour chaque frequence /, & la multiplication de la fonction de directi- 
vite F*\0 9 0j) par la fonction de directivity G(0 9 0). 

Dans le cas ou le traitement est un traitement de type convolution, les 
constantes d t ^ r%m . sont d6termin6es par la relation : 



J',m' 



= < 



A « SrS* + <-m,w) Pout TU > 0 



J'.m' (1) 



J (1) 
c W',m" 



pour m = 0 



J_ (J** W _ J',m' (1) ^ nmir m < n 
, y/2 V C I,-m,l",m" c J,m,J",m",/ P° ur 771 < U 



avec 



J'.m' (1) 
°J,m,f",m" 



1 / J',m' (2) , /',m' (2) \ // ^ r> 

J\m' (2) w n 

C <,m,/",0 POUr ™ = 0 

J_ f/>™' ( 2 ) _ J'*™' (2) \ nonr _// < A 



avec 



<',m' (2) _ 



1 (I'sn* (3) <',-m' (3H m' > 0 



J'fi (3) 



pour m! = 0 



1 / /',-m' (3) _ *',m' (3) \ , fl 



avec 
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I'.m' (3) _ ri rm'-m o/T^m' r m " V **" 1 '+M ~f+|m'| ^' r* 

\| ^ U t'+|m»| *=|m| 

min(i'-m',J'+m") ( /np+k 

Y: (-ir k+m c?, +m „ ctt x i 1/c ^ 51 m ^ 



p=max(0,-m) * '* l/C?,+? +m si m < 0 

et 

eZ = (signe(m)) T 



5 Ce traitement correspond d I'application d'un filtrage spatial avec fil- 

trage temporel. 

En effet, les coefficients d'harmoniques sph6riques Gi^JJ) corres- 
pondent & des coefficients de Fourier-Bessel d'un champ acoustique produit par 
un systeme acoustique lin£aire et invariant angulairement, en reponse & une ex- 

10 citation par un champ acoustique de reference tel qu'une onde plane impulsion- 
nelle provenant de la direction (0,0). Une telle r6ponse est designee par 
['expression « reponse spatio-temporelle ». Un systeme acoustique lineaire est 
invariant angulairement lorsque qu'une rotation du champ qu'il regoit induit la 
meme rotation sur le champ qu'il g6n6re. Ainsi, les coefficients G lim {j) corres- 

1 5 pondent a la reponse spatio-temporelle d'un tel systeme acoustique. 

Les coefficients de Fourier-Bessel P^iM) du champ acoustique initial 
sont 6galement considers comme des coefficients d'harmoniques spheriques 
d'une fonction de directivity P^{0^J). Ainsi, la convolution spatiale avec filtrage 
temporel se traduit, pour chaque frequence f 9 par une convolution spatiale de la 

20 fonction de directivity l*\0 9 #j) par la fonction de directivity G(0,tff) 9 qui est consi- 
d6r6e & chaque frequence / comme un noyau de convolution directif. 

Dans le mode de realisation dgcrit, la convolution est appliquee en 
consid6rant une direction (#»#)), le noyau de convolution G{0^J) est alors tourn6 
en amenant son pole nord (direction verticale avec 0- 0 et ^= 0) en (#>,#)). Pour 

25 obtenir la valeur du champ acoustique modifte I^IJK+fi dans la direction (^,^0, 
la valeur dans la direction du noyau tourn6 est multipliee avec la valeur de 
F^iQffl dans |a direction (#),{&) et ces produits sont somm6s pour toutes les di- 
rections (#>,#>). 

Ce traitement appliqu6 au champ acoustique permet, par exemple, de 
30 modifier l f 6talement angulaire, soit la pr6cision de localisation spatiale, c'est-a- 
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dire la ponctualite des sources virtuelles, dont on considSre qu'elles generent le 
champ acoustique initial, pour chaque frequence / Ce traitement est I'analogue 
spatial de la convolution temporelle par une fenetre, qui a pour effet de modifier 
T6talement temporel d'un evenement sonore, comme par exemple un effet de 
moyennage ou de retard. 

Dans le cas ou la fonction de directivity G l/n {J) est a symetrie cylindri- 
que, c'est-S-dire si G/,m(/)=0 pour m*0 f les equations ci-dessus se simplifient, et 
les coefficients cfiSv sont alors d6termin6s par la relation : 



,m" = S l" S l S m' 6 m" \j 



47T 

2/ + 1 



La fonction de directivity G(0 9 0j) peut etre ind6pendante de la fre- 
quence et est alors note G(0 9 $. Le traitement est alors uniquement spatial et 
correspond, pour chaque frequence/, a la convolution de la fonction de directivite 
P^\^J) associee au champ initial par la fonction de directivite G(0,0). 

Par ailleurs, le proc6d§ de I'invention permet 6galement la r6alisation 
de traitements de rotation. Dans ce cas cependant, il n'est pas fait usage de fonc- 
tions de directivite et le precede requiert des param§tres (0,0, yd, repr6sentatifs 
d'une rotation du champ acoustique autour de Torigine O. 

Ainsi, si le traitement est une rotation, les combinaisons de filtrage 
sont determin6es d partir des param6tres (0,0,yd representatifs de la rotation du 
champ acoustique selon la relation : 



avec 



U m'm 



°/,m — °l U rn'sn 

$t(D l m ,, m + DL m ,, m ) pourm' 

V^»(lCo) pourm' 

&(/^m + #-m',m) pourm' 

y/2R(D l 0tTn ) pourm' 



pour m 
pour m' 
pour m! 
pour m' 
pour m' 



> 0 et 

> 0 et 

> 0 et 
= 0 et 
= 0 et 
= 0 et 

< 0 et 

< 0 et 

< 0 et 



m > 0 
m = 0 
m < 0 
m > 0 
m = 0 
m < 0 
m > 0 
m = 0 
m < 0 



avec 
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D l m >, m = C C: <4', m (*) e"' m '* 

avec 



//; 7 /, Za] min(Z-m',i+m) 

<*~ w - V (i + «>■ (i - «)■ _ir m - m , ( 5 ,+ra 

/ 0\ 2/+m-m'-2/c ✓ ^ \ 2k+m'—m 

C cos 2 ; ( sin 2j 



et 

5 5R(x) = partie reelle (x) 

et 

S(x) = partie imaginaire (x) 



Dans cette relation la convention utilisee pour les angles 0, 0 et y/> 
10 consiste & effectuer une rotation d'un angle ^autour de I'axe (OZ) puis d'un an- 
gle 0 autour de I'axe (OY) puis d'un angle ^ autour de I'axe (OZ). 

Enfin, le traitement effectu6 peut 6galement correspondre £ un traite- 
ment neutre. Dans ce cas v les combinaisons de filtrage sont toutes nulles, sauf 
pour les valeurs ou /*= / et m'= m ou elles valent 1 . 
15 Bien entendu, d'autres traitements ainsi que des combinaisons de trai- 

tements decrits peuvent etre r6alis6s grace & la definition de I'ensemble de cons- 
tantes cfg^ de pond6ration. 

Par exemple, un traitement a!6atoire peut etre r6alis6 par la definition 
al6atoire de I'ensemble de constantes. 
20 Avantageusement, la sous-6tape 30 de calcul est conduite uniquement 

pour les constantes c$;Z? rM non nulles. 

Les combinaisons de filtrage Cf£'(f) ainsi obtenues sont plac6es dans 

une matrice Cde taille (Z/+1) 2 sur (L+\) 2 . 

Dans le cas ou le traitement a appliquer consiste uniquement en une 
25 rotation, une convolution, ou un filtrage, la matrice C correspond directement aux 
combinaisons de filtrage d appliquer et est notee de mantere g6n6rale Fde sorte 
que les resultats de la sous-etape 30 de calcul forment les combinaisons de fil- 
trage destinies d etre appliqu6es aux coefficients de la representation initiate. 

Dans le cas ou le traitement correspond d une combinaison de traite- 
30 ments tels que des rotations, filtrages, convolutions et/ou traitements quel- 
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conques, la sous-etape 30 est r6p6t6e une pluralite de fois et chaque matrice C 
est successivement memorisee afin d'etre combin6e avec les matrices suivantes 
pour finalement deiivrer une matrice F comprenant les coefficients correspondant 
aux combinaisons de filtrage destinees & etre appliquees aux coefficients de la 
5 representation initiale. Ainsi, la matrice F est obtenue par multiplication matricielle 
des matrices C successives. 

Chaque element F^fif) de la matrice F represente la reponse en fre- 
quence d'un filtre de traitement spatio-temporel. 

Diff6rentes formes de filtres et de methodes de filtrage peuvent etre 
10 mises en oeuvre lors de I'etape 10 duplication des combinaisons de filtrage a la 
representation initiale du champ acoustique. 

Dans le cas ou les filtres sont param6tr6s directement avec des re- 
ponses en frequence, les coefficients Ffjf(J) sont directement appliques. Avanta- 

geusement, I'etape 6 de determination des combinaisons de filtrage comprend 
15 une sous-etape 32 de conversion afin de determiner les parametres des filtres 
pour d'autres methodes de filtrage. 

Par exemple, les combinaisons de filtrage F t % m \f) sont converties en : 

- des reponses impulsionnelles finies calcuiees par transform6e de 

Fourier temporelle inverse de F^{f) % chaque r6ponse impulsionnelle 

20 est 6chantillonnee puis tronqu6e d une longueur propre a chaque 

reponse ; ou 

• des coefficients de filtres recursifs d reponses impulsionnelles infinies cal- 
cuiees d partir des Ffiftf) avec des methodes d'adaptation. 

Les combinaisons de filtrage ainsi definies sont appliquees lors de 
25 I'etape 10 d'application des combinaisons de filtrage laquelle transforme la repre- 
sentation initiale /^/^(Z) en une representation modifiee notee F^tjjj) corres- 
pondant & un champ acoustique modifie, en lui appliquant les filtres de traitement 
representatifs du traitement au moins spatial defini lors de I'etape 2. 

Lors de cette 6tape 10, des coefficients p (T) i^(0 correspondant aux 
30 coefficients de la representation modifi6e P (7) , sont deduits des coefficients 
/>%7i(0 correspondant aux coefficients de la representation initiale P ff) % par 
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I'application des filtres de traitement de r6ponse en frequence Ffi?(f) de la fagon 
suivante : 

ou /*\n(/) est la transform^ de Fourier de p (i) /^,(0 et ou P^i^if) est la transfor- 
5 m6e de Fourier de /> (7) /^,(0, coefficients relatifs a la representation initiate et res- 
pectivement £ la representation modifiee. 

Ainsi, I'etape 10 met en oeuvre les methodes classiques de filtrage 
temporel pour ('application des combinaisons de filtrages delivrees par retape 2 
dans une forme adaptee. 
10 Parexemple: 

- Le filtrage peut etre effectue au moyen de methodes de filtrage dans 
le domaine frequentiel, comme par exemple des techniques de convolution par 
bloc. Dans ce cas les combinaisons de filtrage sont foumies sous la forme des 
r6ponses en frequences Ffiftf) . 

15 - Le filtrage peut etre effectue dans le domaine temporel par convolu- 

tion temporelle. Dans ce cas, les combinaisons de filtrage sont foumies sous la 
forme des reponses impulsionnelles f*£\t) . 

- Le filtrage peut etre realise dans le domaine temporel au moyen de 
filtres r6cursifs & reponse impulsionnelle infinie. Dans ce cas, les combinaisons 

20 de filtrage sont foumies sous la forme des coefficients des relations de recur- 
rence. 

Sur la figure 4, on a d6crit le cas d'un filtrage par reponse impulsion- 
nelle finie. Ce filtrage requiert initialement la determination d'une variable 7£ m ' 

correspondent au nombre d'echantillons propres £ chaque reponse f'jfXQ ■ ce 
25 qui conduit d I'expression de convolution temporelle suivante : 

/=0 m=-/ r=0 

Ces coefficients /? (7) /^(0 sont repr6sentatifs d'un champ acoustique 
correspondant au champ acoustique initial aprds I'application du traitement au 
moins spatial qui associe a un champ acoustique produit par au moins une 
30 source, un autre champ acoustique produit par au moins une source dont la di- 
rection, reioignement et/ou le gain sont ainsi modifies. 
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Dans un autre mode de realisation, le precede de I'invention permet en 
outre 1'application de traitements de type distorsion correspondent a une distor- 
sion de parties du champ acoustique par des modifications des caracteristiques 
d'etalement, de position, d'amplitude de tout ou partie des sources virtuelles ge- 
5 neratrices du champ acoustique initial. 

Ce traitement de distorsion peut varier avec la frequence, ce qui per- 
met d'effectuer des modifications d'6loignement et de spectre de tout ou partie de 
ces sources virtuelles. 

Avantageusement, ce traitement peut etre combine a un des traite- 
10 ments tels que d6crits en reference a la figure 3. 

Sur la figure 5, on a d6crit Torganigramme du proced6 de I'invention 
dans un tel mode de realisation combinant un premier traitement appeie traite- 
ment nominal, correspondant £ un traitement tel que d6crit pr6c6demment, avec 
un traitement de distorsion spatiale qui se traduit par des distorsions de caracte- 
15 ristiques de frequence, d'amplitude et/ou de phase de composants de la repre- 
sentation du champ acoustique initial. 

Ainsi, dans ce mode de realisation, I'etape 2 de definition du traitement 
comprend, de meme que prec6demment, retape 4 d'etablissement d'un ensem- 
ble de fonctions de directivite comportant la sous-6tape 12 de determination di- 
20 recte et/ou les sous-etapes 14 de saisie et 16 d 'interpretation ainsi que la sous- 
etape 18 de conversion permettant de deiivrer un ensemble de fonctions de di- 
rectivite sous la forme d'un ensemble de coefficients correspondents a la decom- 
position de ces fonctions sur une base d'harmoniques sph6riques. 

Dans ce mode de realisation, une partie de I'ensemble de fonctions de 
25 directivite determine lors de retape 4 correspond d N couples de distorsion repre- 
sentatifs d'une distorsion a appliquer. De maniere similaire au mode de realisa- 
tion d6crit precedemment, ces couples de distorsion sont d6crits sous la forme 
d'un ensemble de coefficients, notes {(Hi^Jfijr Wtf))}. correspondant & leur 
decomposition sur une base d'harmoniques spheriques. 
30 Chaque couple de distorsion est constitue de deux fonctions de direc- 

tivite qui correspondent aux coefficients de Fourier-Bessel d'un champ acousti- 
que source et d'un champ acoustique cible, l'op6ration de distorsion applique au 
champ acoustique source tendant d le faire ressembler au champ acoustique 
cible. 
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Ainsi, retape 4 delivre un ensemble de fonctions de directivite corres- 
pondant au traitement nominal ainsi qu'e des couples de distorsion representatifs 
de la distorsion £ appliquer sous la forme d'un ensemble de coefficients corres- 
pondant a la decomposition des fonctions de directivite sur une base 
5 d'harmoniques sph6riques. 

De m§me que pr§c6demment, retape 2 comporte egalement retape 
20 de determination de I'ensemble de parametres du traitement par le biais de la 
sous-etape 22 de determination directe et/ou des sous-etapes 24 de saisie et 26 
d'interpretation. 

10 Toutefois, dans ce mode de realisation, I'etape 20 de determination de 

parametres permet en outre la determination d'un parametre A, pouvant §tre 
ajuste entre 0 et 1, specifiant la ressemblance entre le traitement final a appliquer 
et le traitement nominal determine par retape 4, soit encore rimportance relative 
de la distorsion a effectuer. 

15 Ainsi, si A est 6gal £ 1, seule la distorsion est prise en compte et si A 

est egal & 0 f seul le traitement determine lors de retape 4 est pris en compte. 

De meme que dans le mode de realisation decrit en reference & la fi- 
gure 3, retape 6 de determination des combinaisons de filtrage comprend la 
sous-etape 30 de calcul realis6e ainsi que decrit precedemment et permettant de 

20 deiivrer une matrice C comprenant un ensemble de combinaisons de filtrage 
correspondant d ('application du traitement nominal tel que defini par I'ensemble 
de fonctions de directivite determine & retape 4. 

L'ensemble de couples de distorsion est utilise dans Tetape 6 de de- 
termination des combinaisons de filtrage lors d'une sous-etape 50 de distorsion 

25 par modification des combinaisons de filtrage d6termin6es £ Tissue de la sous- 
etape 30 de calcul. 

Lors de cette sous-etape 50, les combinaisons de filtrage Cf£'(f) , sont 
modifi6es d partir des couples de distorsion notes {{HijnJfiJT i^M))} e * du para- 
metre A afin de foumir en sortie les combinaisons de filtrage Fffif) representati- 

30 ves du traitement incluant le traitement nominal et la distorsion et destin6es § 
etre appliqu6es aux coefficients de la representation initiate. 

Lors de la sous-etape 50, les coefficients d'harmoniques spheriques 
Hijnjj) representatifs des N fonctions de directivite sources sont places dans une 
matrice Hde taille (L+l) 2 sur Nde la fagon suivante : 
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#o,o,i(/) H 0.0.2(f) ■ 
.-,.,(/) #,.-,. 2 (/)- 
#>,<>..(/) #..o. 2 (/)- 
#>....(/) H UU2 (f)- 



H L 'o,M H L ^. 2 (f) - 
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' Hq.q.nV) 
#..-..*(/) 
■ #..<>.*(/) 



De meme, les coefficients d'harmoniques spheriques /T/^(/) repre- 
sentatifs des N fonctions de directivite cibles sont places dans une matrice IT de 
taille (L'+l) 2 sur N : 



H'o.o.iif) H'o.oaif) • 

#'..-.,.(/) H\,. it2 (f)- 

**\,o.i(f) H\.o.2(/)- 

#',....(/) H\ A M- 



H' v ,. L , x (f)H\ % . va {fy- 
H\ •.?,•(/) #'r.p.2(/) " 



'H\\-l\n(J) 
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Les coefficients Cf£'(f) des combinaisons de filtrage du traitement no- 
minal sont egalement places dans une matrice Cde taille (L'+l) 2 sur (Z.+1) 2 : 



CS:SV) c^V) c,V(/) c?fU) - 
cis'c/) c.'.-.-.'c/) c,V(/) c.vr'c/) - 
ci:S(/) c.'/f.c/) c,Vo°(/) c,',f(/)- 
c^c/) c.'.J.c/) c,v P '(/) c.'.vc/) .. 


Ci;;I(/)... 


C2$(f) ••• 
Cis«(/)- 
Ci;S(/)- 
Ci;i(/)-- 


ci-XU) 
civvy 
ci%V) 
ciiV) 




ci'#'if>>- 




ct:i L \f) 


c 0 £ ,i°(/) c.^c/) c,^°c/) c,^ 0 (/) .. 


ct'^if)- 


c t$(f>- 


ct'fV) 


.Co^x/) c.^Vc/) c.yx/) cwxn- 




C &'(/)- 


ctixj) . 



Les combinaisons de filtrage /#T'(/) representatives du traitement 

dans son ensemble sont ensuite d6termin6es au cours de la sous-etape 50 de 
15 distorsion correspondant d une optimisation qui determine un ensemble de com- 
binaisons de filtrage v6rifiant au mieux les couples de distorsion et la ressem- 
blance entre F t % m \f) et C^ij) . 
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Cette optimisation peut etre r6alis6e au moyen de methodes telles que 
des methodes de descente de gradient ou diversion matricielle et permet ega- 
lement de determiner, sur les sources virtuelles dont on considere qu'elles gene- 
rent le champ acoustique initial, des modifications qui ne sont pas directement 
5 specifies par les couples de distorsion. 

Par exemple, les modifications d6finies par les couples de distorsion 
sont interpolees sur les parties voisines des parties directement visees de ma- 
niere a obtenir les combinaisons de filtrage, ce qui se traduit par des modifica- 
tions progressives de la representation initiale du champ acoustique autour des 
10 parties directement visees par les couples de distorsion. 

Le mode de realisation decrit utilise une m6thode d'inversion matri- 
cielle qui conduit & la determination d'une matrice F lors de la sous-etape 50 se- 
lon la relation : 

F = ( (l-A) C+JUTlf) ( (1-^) / + AHH T ) 1 
15 ou / designe la matrice identity de taille (Z+l) 2 . 

Dans le cas ou X vaut 1 et ou N< (Z+l) 2 , la matrice F est donn6e par : 
F= C + (IT - CH) ( H T H) 1 H T 



Dans la matrice F, les coefficients Ffj?(f) sont organises de la ma- 
20 niere suivante : 



F 0 °6°(f) F«*(f) Ftfif) F«/(f) F L °;° L (/)... F L ° L Q (f) 

F^if) F^if) F t \f '(/)"• F L \:l<f) -- Fl;o l (f) ~ F^(J) 

FSAif) W(/) Ftfif)- F£;« L (f) ~ F£$(f) ■» Ftfif) 

F 0 \b!(f) F^\(J) FtfV) F^{f) Fl;l L {f) ~ F^(f) F L \l(f) 

• • « . • . . 

FoWif) F£f(f) F^fif) F™(f) - F^if)" Ftfif)- Ftfif) 

• • • » . • • 

Fo&Xf) Ftf Xf) F X &V) F^V)--- F^if '(/)... F&y) ~ F&Xf) 



Bien entendu, selon la methode d'interpolation utilis6e f cette etape 
25 d'optimisation peut entramer des variations progressives ou au contraire brutales 
des modifications entre les parties directement visees par les couples de distor- 
sion et les parties voisines. 
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Par la suite, les combinaisons de filtrage sont 6ventuellement 

converties lors de la sous-etape 32 et sont appliquees lors de Tetape 10 ainsi que 
cela a 6t6 decrit en reference aux figures 3 et 4. 

Avantageusement, plusieurs distorsions sont appliqu6es successive- 
5 ment, le proc6d6 comportant alors une plurality de sous-etapes 50 de distorsion, 
recursives les unes par rapport aux autres, de sorte que chaque sous-etape 50 
regoit en entree, comme traitement nominal, les combinaisons de filtrage deii- 
vrees par la sous-etape pr6c6dente. 

Dans le cas ou le parametre X vaut 0, la distorsion n'est pas appliquee 
10 et les combinaisons des filtrages dgterminees lors de la sous-etape de calcul 30 
correspondent directement aux combinaisons de filtrage representatives du trai- 
tement. Ainsi parametre, ce mode de realisation correspond au premier mode de 
realisation d6crit en reference aux figures 3 et 4. 

Dans un autre mode de realisation, seule la distorsion est appliqu6e et 
15 les combinaisons de filtrages sont determin6es uniquement au moyen des cou- 
ples de distorsion selon la relation : 
-SiAT<(L+l) 2 : 

F = I + (IT - H) ( H T H ) A H J 

-SiW>(Z,+i) 2 : 
20 F = IT H r ( hh t y 1 

Dans ces equations, les matrices H et IT sont determinees ainsi que 
cela a ete decrit precedemment. Dans ce cas, la matrice C n'est pas utilis6e et la 
sous-etape 30 n'est pas r6alis6e. Ainsi, retape 2 comporte uniquement la sous- 
etape 50 de distorsion. 
25 Ce mode de realisation est une simplification du mode de realisation 

decrit lorsque X vaut 1 et que le traitement nominal est le traitement neutre. 

II apparaTt done que le procede de I'invention permet d'appliquer des 
traitements quelconques, pr6sentant une tres large diversite d la representation 
d'un champ acoustique presentant un haut niveau de precision grdce a 
30 Tutilisation d'un ensemble de fonctions de directivite representatif du traitement et 
repr6sent6 sous la forme d'un ensemble de coefficients correspondant d la de- 
composition de ces fonctions sur une base d'harmoniques spheriques. 

Dans le cas ou le traitement d appliquer peut se decomposer en une 
plurality de traitements 6l6mentaires tels que d6crits precedemment, les combi- 
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naisons de filtrage du traitement dans son ensemble peuvent etre obtenues par 
multiplication matricielle des matrices correspondant d chacun de la plurality de 
traitements, soit encore par arrangement des combinaisons de filtrage d6livr6es 
par la sous-6tape 30 de calcul et/ou la sous-etape 50 de distorsion. 
5 Ce traitement est ainsi r6pete une pluralite de fois, les r6sultats de 

chaque occurrence etant combines avec ceux de ('occurrence precedente. 

Sur la figure 6, on a represents un schema synoptique d'un dispositif 
adapts pour la mise en oeuvre du procede de invention. 

Ce dispositif 60 permet d'effectuer un traitement au moins spatial 
10 d'une representation d'un champ acoustique afin de deiivrer une representation 
modifiee par le traitement au moins spatial. 

Le dispositif 60 comporte des moyens 61 de definition du traitement 
pour la mise en oeuvre de I'etape 2 du procede de Tinvention. Ces moyens 61 
comportent des moyens d'etablissement d'un ensemble de fonctions de directivi- 
15 te repr6sentatif du traitement, sous la forme d'un ensemble de coefficients 
correspondant & la decomposition de ces fonctions de directivite sur une base 
d'harmoniques spheriques. 

Avantageusement, ils comportent un module 62 de saisie directe, par 
un op6rateur, & partir d'une memoire ou d'un dispositif exterieur, d'informations 
20 de fonctionnement du dispositif 60 permettant notamment la determination di- 
recte de I'ensemble de fonctions de directivite. 

Par exemple, ce module 62 de saisie est constitue d'un clavier permet- 
tant la saisie directe d'informations ou encore d'une interface avec un dispositif 
exterieur fournissant directement les informations n6cessaires d la mise en oeu- 
25 vre du procede de I'invention. 

Le module 62 peut 6galement etre mis en oeuvre pour Tentr6e d'une 
consigne de traitement, les moyens 61 comportent alors egalement un module 64 
d'interpretation automatique de cette consigne, afin de deiivrer des informations 
de fonctionnement du dispositif 60 et notamment I'ensemble de fonctions de di- 
30 rectivite. 

Les informations entrees grace au module 62 de saisie et/ou 64 
d'interpretation, peuvent en outre etre de diff6rents formats predetermines et le 
dispositif 60 comporte un module 66 de conversion des informations entre des 
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formats predetermines et des formats selectionn6s pour la mise en oeuvre du 
proc6d6 de I'invention. 

Notamment, le module 66 permet de convertir ledit ensemble de fonc- 
tions de directivity de Tun ou I'autre d'une plurality de formats en un ensemble de 
5 coefficients correspondant & leur decomposition sur une base d'harmoniques 
sph&riques. 

Ainsi, les modules 62, 64 et 66 permettent la mise en oeuvre des sta- 
pes 4 et 20 du proc6d6 de I'invention, de manure a delivrer un ensemble de 
fonctions de directivity representatif du traitement sous la forme des coefficients 
10 correspondant a la decomposition de ces fonctions sur une base d'harmoniques 
sph6riques. 

Dans le cas ou le dispositif 60 est utilise pour la realisation d'un traite- 
ment de distorsion, les modules 62, 64 et 66 deiivrent des fonctions de directivite 
formant des couples de distorsion. 
1 5 Les modules 62, 64 et 66 permettent egalement la determination de 

tout ou partie d'un ensemble de parametres necessaire a la realisation du proce- 
d6 de Tinvention. 

Notamment, ces modules 62, 64 et 66 permettent la definition de tout 
ou partie des parametre suivants : 
20 - le parametre L sp6cifiant I'ordre limite de fonctionnement du proc6de 

et correspondant d la limite mathematique de precision souhaitee ; 

- un parametre S definissant la nature du traitement £ effectuer ; 

- tout ou partie de I'ensemble de constantes not6es c££.,m- formant 

des coefficients de ponderation et definies en fonction de la nature du traitement ; 
25 et 

- le parametre X % sp6cifiant I'importance d'une distorsion a effectuer. 
Avantageusement, le dispositif 60 comporte une memoire permettant 

de stocker tout ou partie de ces parametres et notamment I'ensemble de cons- 
tantes. Par exemple, une telle memoire comporte une pluralite d'ensembles de 
30 constantes, chacun correspondant £ un traitement, I'ensemble de constantes mis 
en oeuvre etant selectionne dans cette pluralite en fonction du parametre S re- 
presentatif du traitement. 
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Alternativement, les moyens 61 de definition du traitement compren- 
nent un module de calcul adapts pour determiner tout ou partie de I'ensemble de 
constantes en fonction du parametre S representatif de la nature du traitement. 

L'ensemble de ces informations de fonctionnement est introduit dans 
5 des moyens 70 de determination de combinaisons de filtrage, compris dans les 
moyens 61 de definition du traitement, et adaptes pour la mise en oeuvre de 
I'etape 6 du precede. Ces moyens 70 comportent un module 72 de calcul par 
application de combinaisons Iin6aires permettant la mise en oeuvre de la sous- 
etape 30 de calcul du procedg. Ce module 72 regoit en entree I'ensemble de 
10 fonctions de directivite afin de deiivrer des combinaisons de filtrage par combi- 
naisons Iin6aires des coefficients correspondants a la decomposition sur une 
base d'harmoniques spheriques des ensembles de fonctions de directivite repre- 
sentatives du traitement. 

Avantageusement, les constantes deiivrees par les modules 62, 64 et 
15 66, ferment des coefficients de pond6ration et sont utilis6es comme coefficients 
des combinaisons Iin6aires dans le module 72. 

Dans le cas ou le dispositif 60 met egalement en oeuvre un traitement 
de distorsion, il comporte egalement un module 74 de distorsion par modification 
des combinaisons de filtrage deiivrees par le module 72, mettant en oeuvre la 
20 sous-etape 50 du proc6d6. Ce module regoit en entr6e les couples de distorsion 
repr6sentatifs de la distorsion a appliquer, ainsi que le parametre X representatif 
de Teffort de distorsion mis en oeuvre ainsi que cela a 6te d6crit pr6cedemment. 

Dans le cas ou seule une distorsion est appliquge, le module 72 
n'intervient pas dans la mise en oeuvre du dispositif et uniquement le module 74 
25 de distorsion est n6cessaire. 

Ainsi, les moyens 70 deiivrent un ensemble de coefficients corres- 
pondant d des combinaisons de filtrage representatives du traitement a appliquer. 

Eventuellement, ces combinaisons de filtrage sont introduites dans un 
module 76 de conversion adapte pour deiivrer des parametres de filtres & partir 
30 des combinaisons de filtrage en mettant en oeuvre la sous-etape 32 du precede 
ainsi que cela est d6crit prec6demment. 

Ces parametres de filtres ou directement les combinaisons de filtrage 
sont appliques £ des moyens de filtrage 80 qui mettent en oeuvre retape 10 du 
procede et regoivent en entr6e la representation initiate du champ acoustique 
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sous forme d'un ensemble de coefficients representatif de ce champ dans le 
temps et dans les trois dimensions de I'espace pour appliquer les filtres determi- 
nes afin de dglivrer une representation modifiee correspondant au champ acous- 
tique initial modifie par le traitement defini pr6c6demment. 
5 Par exemple, les moyens de filtrage 80 sont des moyens de filtrage 

num6riques appliques par un processeur de traitement du signal d une represen- 
tation numerique du champ acoustique initial. 

Dans d'autres modes de realisation, ces moyens de filtrage sont des 
moyens de filtrage analogiques d gain commande, appliques a une repr£senta- 
10 tion analogique du champ acoustique sous la forme d'un ensemble de signaux 
correspondant chacun d un coefficient de la representation initiate du champ 
acoustique. 

Par exemple, le dispositif 60 peut §tre realise sous la forme d'un ordi- 
nateur 6quipe de moyens logiciels de calcul et de moyens materiels d'acquisition 

15 et de restitution, de maniere a permettre la mise en oeuvre des moyens 62 
d'entree d'informations ainsi que ('acquisition de la representation initiate du 
champ acoustique et la sortie de la representation modifiee du champ acousti- 
que, les diff6rents elements du dispositif 60 etant realises par des combinaisons 
de moyens logiciels et de moyens materiels. 

20 Dans d'autres modes de realisation, le dispositif 60 correspond d une 

carte d'acquisition sp£cifique dediee uniquement d la realisation de traitements 
acoustiques ou encore sous la forme d'un composant programme pour la mise 
en ceuvre de ces fonctions uniquement. 

Avantageusement, le dispositif 60 est tnt6gr6 dans un systeme acous- 

25 tique comprenant des moyens d'acquisition d'un champ acoustique et des 
moyens de representation de ce champ acoustique sous la forme d'un ensemble 
de coefficients dans le temps et dans les trois dimensions de I'espace afin de 
deiivrer la representation initiale du champ acoustique. Un tel systeme comporte 
6galement un ensemble de restitution commande par des moyens de pilotage 

30 recevant en entree, la representation modifiee du champ acoustique et pilotant 
un ensemble de restitution afin d'obtenir la restitution de ce champ acoustique 
modifie. 
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REVENP1CATIONS 

1 . Precede de traitement d'une representation d'un champ acoustique, 
afin de deiivrer une representation modiftee par application de ce traitement 
correspondant a un traitement au moins spatial du champ acoustique, compor- 

5 tant : 

- retablissement (8) d'une representation initiate du champ acoustique 
sous la forme d'un ensemble de coefficients repr6sentatifs de ce champ dans le 
temps et dans les trois dimensions de I'espace ; et 

- retablissement (10) d'une representation modifi6e par des combinai- 
10 sons de filtrage correspondant audit traitement, appliquees auxdits coefficients de 

ladite representation initiale, 

caracterise en ce qu'il comporte en outre une etape (2) de definition 
dudit traitement, laquelle comprend : 

- une etape (4) d'etablissement d'un ensemble de fonctions de directi- 
15 vite representatif dudit traitement, sous la forme d'un ensemble de coefficients 

correspondant a la decomposition de ces fonctions sur une base d'harmoniques 
spheriques ; et 

- une etape (6) de determination desdites combinaisons de filtrage a 
partir de cette decomposition pour la mise en oeuvre d'une etape (10) 

20 d'application des combinaisons de filtrage ainsi definies auxdits coefficients de la 
representation initiale, afin d'obtenir ladite representation modifiee. 

2. Precede selon la revendication 1, caracterise en ce que ladite etape 
(4) de determination d'un ensemble de fonctions de directivite comporte une 
sous-etape (12) de determination directe de tout ou partie dudit ensemble de 

25 fonctions de directivite dans I'un ou I'autre d'une pluralite de formats predetermi- 
nes. 

3. Procede selon I'une quelconque des revendications 1 ou 2, caracte- 
rise en ce que ladite etape (4) de determination d'un ensemble de fonctions de 
directivite, comporte une sous-etape (16) d'interpretation automatique d'une 

30 consigne de traitement permettant de deiivrer automatiquement tout ou partie 
dudit ensemble de fonctions de directivite dans I'un ou I'autre d'une pluralite de 
formats predetermines, d partir de ladite consigne de traitement. 

4. Proc6de selon I'une quelconque des revendications 1 & 3, caracteri- 
se en ce que ledit ensemble de fonctions de directivite representatif du traitement 
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est exprime selon Tun ou I'autre d'une plurality de formats predetermines et en ce 
que ladite 6tape (4) d'etablissement d'un ensemble de fonctions de directivity 
comporte une sous-6tape (18) de conversion, d partir dudit format predetermine, 
en un ensemble de coefficients correspondant a la decomposition desdites fonc- 
5 tions de directivity sur une base d'harmoniques spheriques. 

5. Proc6d6 selon Tune quelconque des revendications 1 & 4, caracteri- 
s£ en ce que ladite etape de definition (2) du traitement comprend en outre, une 
etape (20) d'etablissement d'un ensemble de parametres dudit traitement. 

6. Precede selon la revendication 5, caract6ris6 en ce que led it en- 
10 semble de parametres comprend tout ou partie des parametres suivants : 

- un ensemble de constantes (c££ r ,*») formant des coefficients de pon- 
deration et defini en fonction de la nature dudit traitement ; et 

- un parametre (L) sp6cifiant I'ordre limite de fonctionnement du pro- 
cede et correspondant a la limite mathematique de precision souhaitee. 

15 7. Proc6d6 selon Tune quelconque des revendications 5 ou 6, caracte- 

ris6 en ce que ladite etape (20) d'etablissement d'un ensemble de parametres 
comporte une sous-6tape (22) de determination directe de tout ou partie dudit 
ensemble de parametres. 

8. Proc6de selon Tune quelconque des revendications 5 e 7, caracteri- 
20 se en ce que ladite etape (20) d'etablissement d'un ensemble de parametres 

comporte une sous-6tape (26) d'interpretation automatique d'une consigne de 
traitement permettant de delivrer automatiquement tout ou partie dudit ensemble 
de parametres d partir de ladite consigne de traitement. 

9. Procede selon Tune quelconque des revendications 6 £ 8, caracteri- 
25 s§ en ce que ladite etape (20) d'etablissement d'un ensemble de parametres 

comprend une sous-etape de calcul de tout ou partie dudit ensemble de constan- 
tes en fonction de la nature dudit traitement. 

10. Proc6d6 selon Tune quelconque des revendications 13 9, caracte- 
rise en ce que ladite etape (6) de determination desdites combinaisons de filtrage 

30 comprend au moins une sous-etape (30) de calcul de combinaisons de filtrage 
par des combinaisons lineaires desdits coefficients correspondant & la decompo- 
sition sur une base d'harmoniques spheriques dudit ensemble de fonctions de 
directivite repr6sentatif dudit traitement. 
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1 1 . Precede selon les revendications 6 et 10 prises ensemble, caracte- 
ris6 en ce que ladite sous-6tape (30) de calcul par des combinaisons lineaires fait 
intervenir ledit ensemble de constantes en tant que coefficients de pond6ration. 

12. Proc6d6 selon Tune quelconque des revendications 10 ou 11, ca- 
5 racteris6 en ce que le resultat de ladite sous-6tape (30) de calcul forme directe- 

ment lesdites combinaisons de filtrage destinees a etre appliquees auxdits coeffi- 
cients de la representation initiate. 

13. Proc6d6 selon la revendication 12, caracterise en ce que ledit trai- 
tement est un traitement de convolution. 

10 14. Proc6de selon la revendication 12, caracterise en ce que ledit trai- 

tement est un traitement de fenetrage. 

15. Procede selon Tune quelconque des revendications 1 d 11, carac- 
terise en ce que ledit traitement comprend au moins une distorsion £ appliquer £ 
la representation initiate du champ acoustique et en ce que tout ou partie des 

15 fonctions de directivite determin6es lors de Tetape (4) d'etablissement d'un en- 
semble de fonctions de directivite ferment un ensemble de couples de distorsion 
repr6sentatif de ladite distorsion sous la forme d'un ensemble de coefficients 
correspondant d la decomposition desdits couples de distorsion sur une base 
d'harmoniques sph6riques. 

20 16. Proc6d6 selon les revendications 10 et 15 prises ensemble, carac- 

terise en ce que ledit traitement comprend un premier traitement dit nominal et la 
ou lesdites distorsions, ladite sous-6tape (30) de calcul deiivrant des combinai- 
sons de filtrage correspondant audit traitement nominal et ladite 6tape (6) de de- 
termination de combinaisons de filtrage comprenant en outre au moins une sous- 

25 etape (50) de distorsion par modification des combinaisons de filtrage deiivrees 
par ladite sous-etape (30) de calcul, a partir desdits couples de distorsion repre- 
sentatifs de la ou desdites distorsions, afin de deiivrer lesdites combinaisons de 
filtrage destinees d etre appliqu6es auxdits coefficients de la representation ini- 
tiate. 

30 17. Proc6d6 selon les revendications 5 et 16 prises ensemble, caracte- 

rise en ce que ledit ensemble de parametres comporte un parametre (X) repre- 
sentatif de I'effort de distorsion, soit encore specifiant la ressemblance entre ledit 
traitement & appliquer et ledit traitement nominal. 
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18. Precede selon la revendication 15, caract6rise en ce que ledit trai- 
tement comprend uniquement la ou lesdites distorsions, ladite etape (6) de de- 
termination des combinaisons de filtrage comprenant au moins une sous-6tape 
(50) de distorsion delivrant des combinaisons de filtrage uniquement £ partir des- 

5 dits couples de distorsion representatifs de la ou desdites distorsions, afin de de- 
livrer lesdites combinaisons de filtrage destinies d etre appliquees auxdits coeffi- 
cients de la representation initiale. 

19. Precede selon Tune quelconque des revendications 16 a 18, ca- 
racterise en ce qu'il comporte une plurality de sous-stapes (50) de distorsion, 

10 recursives les unes par rapport aux autres, de sorte que chaque sous-6tape (50) 
de distorsion regoit en entree les combinaisons de filtrage delivrees par la sous- 
etape prec6dente. 

20. Precede selon Tune quelconque des revendications 16 d 19, ca- 
racteris6 en ce que ladite au moins une sous-6tape (50) de distorsion , permet de 

15 determiner les combinaisons de filtrage correspondant a des modifications & ef- 
fectuer sur des parties de la representation initiale qui ne sont pas directement 
vis6es par lesdits couples de distorsion. 

21. Precede selon Tune quelconque des revendications 10, 16 ou 18, 
caracterise en ce que ledit traitement se decompose en une plurality de traite- 

20 ments elementaires et en ce qu'au moins ladite etape (6) de determination de 
combinaisons de filtrage est rep6t6e une plurality de fois, les combinaisons de 
filtrage resultant des diff6rentes sous-6tapes (30) de calcul et/ou (50) de distor- 
sion et correspondant auxdits traitements elementaires etant arrangees entre 
elles pour former les combinaisons de filtrage destinees a etre appliquees auxdits 

25 coefficients de la representation initiale. 

22. Proc6d6 selon Tune quelconque des revendications 13 21, carac- 
terise en ce que ladite etape (6) de determination des combinaisons de filtrage 
comprend une sous-etape (32) de conversion desdites combinaisons de filtrage 
representees dans un format predetermine pour d6livrer des parametres de filtres 

30 adaptes pour le parametrage de filtres donn6s en vue de la realisation de ladite 
etape (10) duplication des combinaisons de filtrage. 

23. Procede selon Tune quelconque des revendications 1 d 22, carac- 
terise en ce que ledit traitement est un traitement spatial et en ce que ladite etape 
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(10) duplications des combinaisons de filtrage comprend des modifications de 
gain appliquees auxdits coefficients de la representation initiale. 

24. Proc6de selon Tune quelconque des revendications 1 a 23, carac- 
terise en ce que ledit traitement est un traitement spatial et frequentiel et en ce 

5 que ladite 6tape (10) d'application des combinaisons de filtrage comprend 
1'applicatton auxdits coefficients de la representation initiale, de filtres variant en 
fonction de la frequence. 

25. Proc6d6 selon Tune quelconque des revendications 1 a 24, carac- 
terise en ce que lesdits ensembles de coefficients representatifs de la repr6sen- 

10 tation initiale du champ acoustique et de la representation modifiee sont des en- 
sembles de coefficients dits de Fourier-Bessel. 

26. Dispositif de traitement d'une representation d'un champ acousti- 
que afin de deiivrer uhe representation de ce champ acoustique par application 
de ce traitement correspondant £ un traitement au moins spatial du champ 

15 acoustique, ledit dispositif (60) recevant en entr&e une representation initiale du 
champ acoustique sous la forme d'un ensemble de coefficients representatif de 
ce champ dans le temps et dans les trois dimensions de I'espace et delivrant en 
sortie une representation modifiee par les combinaisons de filtrage correspondant 
audit traitement applique auxdits coefficients de la representation initiale, caracte- 

20 rise en ce qu'il comporte des moyens (61) de definition dudit traitement, compre- 
nant : 

- des moyens (62, 64, 66) d'etablissement d'un ensemble de fonc- 
tions de directivite representatif dudit traitement, sous la forme d'un ensemble de 
coefficients correspondant d la decomposition de ces fonctions de directivite sur 

25 une base d'harmoniques spheriques ; et 

- des moyens (70) de determination desdites combinaisons de fil- 
trage a partir de cette decomposition pour leur mise en ceuvre par des moyens 
(80) de filtrage recevant en entr6e lesdits coefficients de la representation initiale 
et delivrant en sortie ladite representation modifiee. 

30 27. Dispositif selon la revendication 26, caracterise en ce que lesdits 

moyens (61) de definition dudit traitement comportent un module (62) de saisie 
directe dlnformations de fonctionnement du dispositif dans Tun ou Tautre d'une 
pluralite de formats. 
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28. Dispositif selon Tune quelconque des revendications 26 ou 27, ca- 
racterise en ce que lesdits moyens (61) de definition dudit traitement comportent 
un module (64) d'interpretation automatique d'une consigne de traitement per- 
mettant de delivrer automatiquement des informations de fonctionnement dans 

5 Tun ou I'autre d'une pluralite de formats predetermines. 

29. Dispositif selon Tune quelconque des revendications 27 ou 28, ca- 
racterise en ce que lesdites informations de fonctionnement comprennent ledit 
ensemble de fonctions de directivite representatif dudit traitement. 

30. Dispositif selon la revendication 29, caracterise en ce que lesdits 
10 moyens (61) de definition dudit traitement comportent un module (66) de conver- 
sion adapte pour convertir ledit ensemble de fonctions de directivite de Tun ou 
Tautre d'une pluralite de formats predetermines en un ensemble de coefficients 
correspondant a la decomposition desdites fonctions de directivite sur une base 
d'harmoniques spheriques. 

15 31. Dispositif selon Tune quelconque des revendications 27 ou 28, ca- 

racteris6 en ce que lesdites informations de fonctionnement comprennent un en- 
semble de parametres dudit traitement. 

32. Dispositif selon la revendication 31 , caracterise en ce que lesdits 
moyens (61) de definition dudit traitement comportent un module de calcul de 

20 tout ou partie dudit ensemble de parametres en fonction de la nature dudit traite- 
ment. 

33. Dispositif selon Tune quelconque des revendications 26 a 32, ca- 
racterise en ce que lesdits moyens (70) de determination desdites combinaisons 
de filtrage comprennent un module (72) de calcul de combinaisons de filtrage par 

25 combinaisons Iin6aires desdits coefficients correspondant & la decomposition sur 
une base d'harmoniques spheriques dudit ensemble de directivite representatif 
dudit traitement. 

34. Dispositif selon les revendications 31 et 33 prises ensemble, ca- 
racterise en ce que ledit ensemble des parametres dudit traitement comprend un 

30 ensemble de constantes et en ce que ledit module (72) de calcul par combinai- 
sons lineaires fait intervenir ledit ensemble de constantes en tant que coefficients 
de ponderation. 

35. Dispositif selon Tune quelconque des revendications 26 d 33, 
caracterise en ce que ledit traitement comprend au moins une distorsion, tout ou 
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partie dudit ensemble de fonctions de directivity formant un ensemble de couples 
de distorsion correspondant a la ou aux distorsions a appliquer d la representa- 
tion initiale et en ce que lesdits moyens (70) de determination desdites combinai- 
sons de filtrage comprennent en outre un module (74) de distorsion par modifica- 
5 tion des combinaisons de filtrage delivrees par ledit module (72) de calcul, a par- 
tir desdits couples de distorsion afin de deiivrer lesdites combinaisons de filtrage 
destinies £ etre appliqu6es auxdits coefficients de la representation initiale. 



racterise en ce que ledit traitement comprend uniquement une ou plusieurs dis- 
10 torsions, tout ou partie dudit ensemble de fonctions de directivity formant un en- 
semble de couples de distorsion correspondant a la ou aux distorsions a appli- 
quer d la representation initiale et en ce que lesdits moyens (70) de determination 
desdites combinaisons de filtrage comprennent uniquement un module (74) de 
distorsion pour deiivrer lesdites combinaisons de filtrage uniquement a partir 
1 5 desdits couples de distorsion. 



36. Dispositif selon Tune quelconque des revendications 26 d 32, ca- 
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